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1.  Introduction 


The  Alekseevski-Tate  equations  (7)  have  long  been  employed  to  predict  the  penetration  of 
kinetic  energy  penetrators  or  rod-like  projectiles.  Typically,  these  equations  are  solved  by 
straight-forward  numerical  integration  techniques  or  by  using  the  exact  solution  developed  by 
Walters  and  Segletes  (2)  and  refined  by  Segletes  and  Walters  (3).  However,  due  to  the  nonlinear 
nature  of  these  equations,  the  exact  solution  yields  penetration  only  as  an  implicit  function  of 
time.  It  is  desirable  to  obtain  an  exact,  albeit  approximate,  solution  to  these  equations  which 
would  yield  the  pertinent  variables,  penetration  velocity,  penetrator  length,  penetrator  velocity, 
and  penetration  as  an  explicit  function  of  time.  The  analytical  nature  of  the  solutions  would 
clearly  reveal  the  interplay  of  the  various  terms  in  the  governing  equations  on  the  solutions. 
Toward  this  end,  as  per  references  (1-2),  a  perturbation  solution  was  obtained  for  the  following 
equation  set: 


if  +  YP  =  (pT  /2)u  2  +  Rt 

(1) 

dv  -  Yp 

(2) 

dt  lpp 

dl 

—  —  u  —v 
dt 

(3) 

dp 

—  or  p  =  \  udt  . 
dt 

(4) 

In  these  equations,  v  is  the  penetrator  velocity,  u  is  the  penetration  velocity,  p  is  the  penetration, 

/  is  the  penetrator  length,  t  is  the  time  after  impact,  Rt  is  the  target  strength  tenn,  YP  is  the 
penetrator  strength  term,  p  represents  the  density,  where  the  subscript  P  stands  for  penetrator  and 
subscript  T  represents  the  target.  First,  the  equations  are  normalized  and  the  method  of 
nonnalization  will  depend  on  the  input  conditions;  namely,  the  p  values,  the  initial  velocity,  and 
the  strength  terms.  For  the  usual  case  of  interest  to  ballisticians  studying  kinetic  energy 
penetrators  impacting  armor  targets,  the  following  normalization  parameters  are  introduced: 

V  =  v/vi ,  U  =  u/vi  ,  X  =  l/L ,  r  =  fit ,  and  (5) 


2 


P 


Pt  /  Pp  •> 


P  = 


V1  +  /V 


(6) 


where  F  is  the  impact  velocity,  L  is  the  initial  penetrator  length,  and  V,  U,  P,  A,  and  r  are  the 
dimensionless  variables,  while  p  ,a,  [1,  and  s  are  the  dimensionless  constants.  The  parameter  ft  is 
used  to  normalize  the  time.  The  constant  s  is  the  perturbation  parameter.  This  is  the  same 


1 


normalization  scheme  used  by  Forrestal  et  al.  ( 4 ).  Forrestal  et  al.  assumed  Rt=0  (i.e.,  a  =  -l) 
for  a  steel  projectile  impacting  a  foundry  core  target  (silica  sand)  at  3.0  km/s. 

Thus,  the  normalized  equations  become 


() V-uf  -2 ccs  =  jU2U2  , 

dV  e  (\  +  ju^ 
dr  ju  ) 


dA 

dr 


'l  +  fP 
v  M  J 


(U-V) 


(7) 

(8) 

(9) 


and 


P  = 


1  +  JU 


T 

J  Udv . 


v  P  j  o 

The  perturbation  method,  as  shown  by  Cole  (5),  involves  letting 


V  =  v0  +eVx  +e2V2  +£3v3  +. 


U  =  U0  +  sUl  +  s2U2  +  £3C/3  +. 


P  —  Pq  +  sP\  +  £  P2  £  Pt,  • 


and  the  perturbation  parameter  £  = 


« 


PpV{ 


1  is  chosen. 


(10) 


(ID 

(12) 

(13) 

(14) 


For  example,  if  a  third-order  perturbation  solution  is  obtained,  terms  of  the  order  of  £to  the 
fourth  power  are  neglected.  Hence,  the  accuracy  of  the  solution  depends  on  the  magnitude  of  s 
and  the  number  of  terms  in  the  above  equation  set. 


One  can  substitute  the  expressions  for  V,  U,  A,  and  P  (equations  1 1-14)  into  the  nondimensional 
equation  set  (equations  7-10)  and  obtain,  to  order  zero  (considering  only  terms  involving  <£°)  ,  the 
following  set: 


{v0-u0y  =m2u20 

n 


(15) 


dV0 

— —  =  0  , 
dr 


(16) 


and 


dA0 

dr 


'ljj? 

V  M  J 


(^0-^0)  . 


where  the  initial  conditions  at  r  =  0  are  Vq  (0)  =  1 ,  A,  (0)  =  1 . 


(17) 
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The  solution  is 


r„=l  , 


A0  =  1  -  r 


po=  —  ’ 


in  agreement  with  Forrestal  et  al.  (4). 

Proceeding  along  these  same  lines  (next  consider  tenns  of  order  e,  etc.),  one  can  calculate  Vo,  V\, 
V2,  V3,  and  the  corresponding  U,  P,  and  X  values.  The  method  is  straight  forward  but  tedious. 

The  final  equations,  to  the  third  order,  are 

V  =  V0+eVl+  e2V2  +  £%  ,  (22) 

U  =  U0  +  sU]  +  e2U2  +  e3U3  ,  (23) 

A  =  Aq  +  +  £  A2  +  £  A3  ,  (24) 


where 


P  —  -Pq  +  «£P|  +  £~  P2  (  £  P3  , 


Vx  =  <91n77  , 


V2  =02n  In 77  +  — — 1  -  — <92[ln 77l2 

77  2  L  J 
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a2  ' 

V, 

77 

-  2n{\n  n)2  +  -i  (in  77)3  -  Q  — +  C3 


C/j  =  —  [in  77 -a] 


A  +  M  fj.  2ju 

=  P^-^Wn)\tlJ_nj_aYn_n  Tfa,  T  t  J 

l  +  /y  n  2  \  2  2  2  )  77  2 1  // v  'J 


=  @[z72  +  77  In  77  ] 
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I  3  2  2  I  1 
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0  2/v  2/U2 
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722  77(ln77)3 
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a2  2a  a2 
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The  normalization  scheme  used  is  deemed  appropriate  for  ballistic  applications,  namely  tungsten 
rods  penetrating  steel.  The  third-order  solution  for  U,  V,  P,  and  X,  with  comparison  to  the  exact 
Alekseevski-Tate  equation  solution,  is  shown  in  figures  1-3.  In  this  case,  Vf  =  2  km/s,  L  =  0.5  m, 
pP=  17,600  kg/m3,  YP=  1.0  GPa,  pT=  7,800  kg/m3,  and  RT=  5.5  GPa.  The  perturbation 
parameter  computes  as  e  =  0.0142,  which  is  much  less  than  one  and  terms  with  a  coefficient  o f  i:4 
were  neglected.  The  agreement  is  excellent  up  to  about  495  ps  for  the  velocities  and  nearly 
identical  for  the  rod  length  and  penetration.  The  exact  (Tate)  solution  terminates  at  about 
518  ps.  The  appendix  lists  the  FORTRAN  program  used  to  calculate  the  rod  velocity, 
penetration  velocity,  rod  length,  and  penetration  all  as  a  function  of  time. 


A  few  comments  are  in  order.  All  equations  contain  a  tenn  like  ln(l-r).  Thus,  a  singularity 
L(  1  +  p) 

occurs  at  r  =  1  or  t  =  — - - =  625.5  ps  for  the  case  described.  Thus,  t  must  be  less  than  this 

juV 

value  for  the  perturbation  solution  to  be  evaluated.  As  this  singularity  is  approached,  deviation 
from  the  exact  solution  occurs.  Hence,  one  can  accurately  calculate  the  perforation  of  finite 
thickness  targets  if  the  penetration  time  is  less  than  this  value.  As  an  example,  figures  4-6  plot 
the  rod  velocity,  penetration  velocity,  penetration,  and  rod  length  for  a  tungsten  rod  impacting  a 
finite-thickness  (25  mm)  steel  plate  as  described,  except  the  initial  rod  length  was  39.12  mm. 

The  first-order  and  third-order  perturbation  theories  are  compared  with  the  Tate  solution.  Again, 
the  agreement  is  excellent.  Note  that  the  first-order  perturbation  theory  also  agrees  very  well 
with  the  exact  solution  for  this  case.  Thus,  approximate  formulae  may  be  derived  for  the 
perforation  of  thin  plates.  The  first  order  equations  are 


V  =  1  +  s&  In  n 


(26) 
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Figure  1.  Comparison  between  rod  and  penetration  velocities  for  a  first-order  perturbation 

solution,  a  third-order  perturbation  solution,  and  the  exact  Tate  solution.  This  case  is 
for  a  tungsten  rod  impacting  a  steel  target  at  an  impact  velocity  of  2000  m/s.  The  initial 
rod  length  was  0.500  m. 


U  =  — —  +  —  [in  77 -a]  , 
1  +  yU  M 


(27) 


A  =  I7  +  £0[n(l-n)+I7\nI7]  , 


(28) 


and 


P  =  — —  [(l-77Xa  +  l)+771n77]  . 


(29) 


When  P=H,  the  finite  target  plate  thickness,  equation  29  can  be  solved  for  the  event  duration  in 
terms  of  77  as 


77  =  exp 


-1 

[H  ,  n-A 

An 

U  7  J 

(1-77) 


77 


(«  +  l) 


A  A  S® 

and  A  =  — 


(30) 
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Figure  2.  Penetration  depth  for  a  first-order  perturbation  solution,  a  third-order  perturbation 

solution,  and  the  exact  Tate  solution.  This  case  is  for  a  tungsten  rod  impacting  a  steel 
target  at  an  impact  velocity  of  2000  m/s.  The  initial  rod  length  was  0.500  m. 

Further,  the  exponential  tenn  in  equation  30  may  be  approximated  by  the  first  two  terms  of  its 
series  expansion  or  exp[x]  =  1  +  x.  In  this  case, 

„  £  +  Vb2+4C  non  1  ,  _  H  1  ,  .  nn 

77  = - where  B=2+a- —  and  C  = - + - (a  +  1)  .  (31) 

2  Aju  ALAju 

From  equation  31,  the  time,  t,  was  25.34  ps  and  the  time  calculated  from  the  more  exact 
equation  30  via  standard  iteration  techniques  was  25.12  ps.  The  rod  velocity  was  1948.156  m/s 
according  to  the  approximate  theory  and  1943.000  m/s  from  the  Tate  solution.  The  penetration 
velocity  was  977.524  m/s  vs.  976.500  m/s  from  the  Tate  solution.  The  final  rod  length  was 
14.20  mm  vs.  14.70  mm  from  the  Tate  solution.  This  close  agreement  means  that  the 
perturbation  solution  is  in  excellent  agreement  with  the  Tate  solution  for  small  times,  even  for  a 
first-order  perturbation  solution.  Forrestal  et  al.  ( 4 )  also  showed  good  agreement  with  Tate  for 
short  times. 
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Figure  3.  Rod  length  for  a  first-order  perturbation  solution,  a  third-order  perturbation  solution,  and 
the  exact  Tate  solution.  This  case  is  for  a  tungsten  rod  impacting  a  steel  target  at  an 
impact  velocity  of  2000  m/s.  The  initial  rod  length  was  0.500  m. 


As  previously  mentioned,  the  perturbation  parameter  must  be  much  less  than  one.  The  value  of  s 
used  in  the  above  equations  implies  a  small  penetrator  strength  or  Y  «  ppV2 .  If  this  is  not  the 

case,  for  example,  for  a  soft  target,  i.e.,  small  RT,  the  equations  can  be  normalized  with 


V  =  v/Vf  ,  U  =  u/Vt  ,  A.  =  l/L,  P 


(32) 


2 


A 


Pt/Pp »  P 


f  \ 
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£  = 


PpVi 


2  ’ 


T  =  fit  , 


and  the  equation  set  becomes 


{v-uf  +2 a£  =  p2u2  , 


(33) 


(34) 
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Figure  4.  Comparison  between  rod  and  penetration  velocities  for  a  first-order  perturbation 

solution,  a  third-order  perturbation  solution,  and  the  exact  Tate  solution.  This  case  is 
for  a  tungsten  rod  impacting  a  25-mm-thick  steel  target  at  an  impact  velocity  of 
2000  m/s.  The  initial  rod  length  was  0.03912  m. 


dV  s 

f1  +  /0 

M 

dr  X 

l  P  J 

,  RT  y 

and 


^lJ\±A{u_v)  . 

dr  y  p  ) 


(35) 


(36) 


This  is  very  similar  to  the  previous  set  of  equations  and  thus  follows  the  same  solution  scheme. 

Note  that  the  normalization  scheme  was  initiated  by  dividing  the  first  equation  by  pp  and  defining 
p  to  be  pT/pp.  Alternately,  one  could  divide  the  equation  by  pT  and  define  p  to  be  pP/pT- 

Other  normalization  schemes  are  possible  namely  if  the  nondimensional  variables  are  defined  as 

*  =  l/L,fi  =Pt/Pp,K  =  ^,Z  =  2(Rt~Yp)  ,  (37) 

Pp  pP 
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Figure  5.  Penetration  depth  for  a  first-order  perturbation  solution,  a  third-order  perturbation  solution, 
and  the  exact  Tate  solution.  This  case  is  for  a  tungsten  rod  impacting  a  25-mm-thick  steel 
target  at  an  impact  velocity  of  2000  m/s.  The  initial  rod  length  was  0.03912  m. 

V  =  U  =  u/y[z,  A  =  *  =  t'^/L-P  =  y/L  •  (38) 

Note  that  Z  >  0  is  required.  Based  on  input,  if  this  is  not  the  case,  one  may  refonnulate  Z7as 

2 (Yp  -Rt) 


The  equation  set  becomes 


(u-v)2  =  JUU2  + 1, 


(39) 


d/ 1 
dr 


=  U -V ,  and 


(40) 


MV  _  -K  _  -Yp  dv  _  £ 

dr  z  2(Rt~yp)  dr  x 


(41) 
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Figure  6.  Rod  length  for  a  first-order  perturbation  solution,  a  third-order  perturbation  solution, 
and  the  exact  Tate  solution.  This  case  is  for  a  tungsten  rod  impacting  a  25-mm-thick 
steel  target  at  an  impact  velocity  of  2000  m/s.  The  initial  rod  length  was  0.03912  m. 

_ _ y 

where  s  = - = - - —  is  chosen  to  be  the  perturbation  parameter.  Again,  the  previous  set  of 

£  2  (Rt-Yp) 

equations  is  very  similar  to  the  original  set  of  equations  and  thus  follows  the  same  solution 
scheme.  The  normalization  scheme  or  equation  set  chosen  will  depend  on  the  input  conditions 
(known  initial  values)  and  the  requirement  to  keep  e  «  1 .  Also,  it  is  advantageous  to  make  the 
time  where  the  logarithmic  singularity  occurs  (x  =  1)  as  large  as  possible. 


2.  Parametric  Studies 


Other  rod-target  configurations  were  studied  using  the  third-order  perturbation  theory  in  order  to 
exercise  the  model.  In  figure  7,  the  rod  velocity  and  penetration  velocity  are  plotted  vs.  time  for 
the  Tate  solution,  the  first-order,  and  the  third-order  perturbation  theory.  In  this  case,  the  initial 
rod  velocity  was  3  km/s.  The  agreement  between  the  third-order  perturbation  model  and  the 
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Figure  7.  Comparison  between  rod  and  penetration  velocities  for  a  first-order  perturbation 

solution,  a  third-order  perturbation  solution,  and  the  exact  Tate  solution.  This  case  is 
for  a  tungsten  rod  impacting  a  steel  target  at  an  impact  velocity  of  3000  m/s.  The 
initial  rod  length  was  0.500  m. 

Tate  solution  is  excellent,  due  to  the  perturbation  parameter  s  =  0.0063,  being  small.  The  length 
and  penetration  are  also  in  excellent  agreement,  as  shown  in  figures  8  and  9.  The  input,  namely 
the  strength  terms,  was  obtained  from  Dr.  Steven  Segletes  of  ARL  for  the  tungsten  rod  vs.  the 
RHA  target. 

Next,  figures  10-12  plot  the  velocities,  penetration,  and  rod  length  vs.  time  for  the  case  studied 
by  Forrestal  et  al.  ( 4 ),  namely  a  steel  rod  impacting  a  geological  (silica-sand)  target.  The  first- 
order  solution,  third-order  solution,  and  the  Tate  solution  are  compared.  The  perturbation 
parameter  was  0.0192  and  the  improvement  of  the  third-order  solution  over  the  first-order 
solution  is  shown.  The  Tate  solution  completes  the  penetration  process  at  56  ps.  The  first-order 
solution  agrees,  within  1%,  up  to  29  ps,  which  is  the  Forrestal  et  al.  ( 4 )  solution.  The  third-order 
solution  is  valid,  again  within  1%,  to  35  ps. 
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Figure  8.  Penetration  depth  for  a  first-order  perturbation  solution,  a  third-order  perturbation 

solution,  and  the  exact  Tate  solution.  This  case  is  for  a  tungsten  rod  impacting  a  steel 
target  at  an  impact  velocity  of  3000  m/s.  The  initial  rod  length  was  0.500  m. 

As  mentioned  previously,  the  strength  terms  for  the  steel  targets  were  obtained  from  Dr. 

Segletes.  Input  values  used  by  Edward  Horwath  of  ARL  were  penetrator  strengths  of  3.6  GPa 
for  tungsten,  3.2  GPa  for  depleted  uranium  (DU),  and  the  target  resistance  was  5.1  GPa  for  RHA. 
All  input  values  used  in  this  study  are  recorded  in  table  1 .  These  tungsten  values,  and  the  value 
previously  given  for  the  RHA  target,  are  designated  as  93%  tungsten  (the  only  difference  in  the 
two  tungstens  is  the  value  assigned  to  Yp)  and  the  plots  are  shown  in  figures  13-15  for  an  initial 
rod  velocity  of  2  km/s.  Again,  the  agreement  is  good.  The  Tate  solution  indicates  that  the 
penetration  process  is  complete  at  652  ps  and  the  third-order  solution  is  valid  up  to  598  ps. 
Figures  16-18  present  the  same  plots  for  an  initial  impact  velocity  of  3  km/s.  In  this  case,  again 
because  of  the  lower  value  of  the  perturbation  parameter,  the  agreement  is  excellent.  The  Tate 
solution  terminates  at  428  ps  and  the  third-order  theory  is  valid  up  to  409  ps.  The  perturbation 
parameter  was  0.05 1 1  for  the  2  km/s  case  and  0.0227  for  the  3  km/s  case.  Table  2  lists  the 
perturbation  parameters  for  all  cases.  Finally,  figures  19-21  give  the  plots  for  an  initial  velocity 
of  2  km/s  for  a  DU  rod  impacting  RHA.  The  agreement  is  good,  with  the  Tate  solution 
terminating  at  640.5  ps  and  the  third-order  solution  valid  up  to  61 1  ps.  The  agreement  is  perfect 
for  the  same  rod/target  configuration  at  an  initial  impact  velocity  of  3  km/s,  see  figures  22-24. 
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Figure  9.  Rod  length  for  a  first-order  perturbation  solution,  a  third-order  perturbation  solution, 
and  the  exact  Tate  solution.  This  case  is  for  a  tungsten  rod  impacting  a  steel  target  at 
an  impact  velocity  of  3000  m/s.  The  initial  rod  length  was  0.500  m. 

Here,  the  perturbation  parameter  was  0.0430  for  the  2  km/s  case  and  0.0191  for  the  3  km/s  case. 
All  input  and  perturbation  values  are  summarized  in  tables  1  and  2. 
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Figure  10.  Comparison  between  rod  and  penetration  velocities  for  a  first-order  perturbation 

solution,  a  third-order  perturbation  solution,  and  the  exact  Tate  solution.  This  case  is 
for  a  steel  rod  impacting  a  silica-sand  target  at  an  impact  velocity  of  3000  m/s.  The 
initial  rod  length  was  0.03912  m. 
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F igure  1 1 .  Penetration  depth  for  a  first-order  perturbation  solution,  a  third-order  perturbation 
solution,  and  the  exact  Tate  solution.  This  case  is  for  a  steel  rod  impacting  a  silica- 
sand  target  at  an  impact  velocity  of  3000  m/s.  The  initial  rod  length  was  0.03912  m. 
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Figure  12.  Rod  length  for  a  first-order  perturbation  solution,  a  third-order  perturbation  solution,  and 
the  exact  Tate  solution.  This  case  is  for  a  steel  rod  impacting  a  silica-sand  target  at  an 
impact  velocity  of  3000  m/s.  The  initial  rod  length  was  0.03912  m. 

Table  1 .  Material  input  properties. 


Material 

Pp  3 

(kg/m  ) 

Pt 

(kg/m3) 

Ye 

(GPa) 

rt 

(GPa) 

Tungsten 

17600 

— 

1.000 

— 

93%  tungsten 

17600 

— 

3.600 

— 

Depleted  uranium  (DU) 

18600 

— 

3.200 

— 

Steel  (RHA) 

— 

7800 

— 

5.500 

Steel  (RHA)  for  93%  tungsten 

— 

7800 

— 

5.100 

Forrestal  penetrator  (steel) 

8000 

— 

1.380 

— 

Forrestal  target  (silica-sand) 

— 

1700 

— 

0.000 
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Figure  13.  Comparison  between  rod  and  penetration  velocities  for  a  first-order  perturbation 

solution,  a  third-order  perturbation  solution,  and  the  exact  Tate  solution.  This  case  is 
for  a  93%  tungsten  rod  impacting  a  steel  target  at  an  impact  velocity  of  2000  m/s.  The 
initial  rod  length  was  0.500  m. 
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Figure  14.  Penetration  depth  for  a  first-order  perturbation  solution,  a  third-order  perturbation 

solution,  and  the  exact  Tate  solution.  This  case  is  for  a  93%  tungsten  rod  impacting  a 
steel  target  at  an  impact  velocity  of  2000  m/s.  The  initial  rod  length  was  0.500  m. 


19 


Figure  15.  Rod  length  for  a  first-order  perturbation  solution,  a  third-order  perturbation 
solution,  and  the  exact  Tate  solution.  This  case  is  for  a  93%  tungsten  rod 
impacting  a  steel  target  at  an  impact  velocity  of  2000  m/s.  The  initial  rod  length 
was  0.500  m. 
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Figure  16.  Comparison  between  rod  and  penetration  velocities  for  a  first-order  perturbation 

solution,  a  third-order  perturbation  solution,  and  the  exact  Tate  solution.  This  case  is 
for  a  93%  tungsten  rod  impacting  a  steel  target  at  an  impact  velocity  of  3000  m/s.  The 
initial  rod  length  was  0.500  m. 
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Figure  17.  Penetration  depth  for  a  first-order  perturbation  solution,  a  third-order  perturbation 
solution,  and  the  exact  Tate  solution.  This  case  is  for  a  93%  tungsten  rod 
impacting  a  steel  target  at  an  impact  velocity  of  3000  m/s.  The  initial  rod  length 
was  0.500  m. 
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Figure  18.  Rod  length  for  a  first-order  perturbation  solution,  a  third-order  perturbation  solution, 
and  the  exact  Tate  solution.  This  case  is  for  a  93%  tungsten  rod  impacting  a  steel 
target  at  an  impact  velocity  of  3000  m/s.  The  initial  rod  length  was  0.500  m. 


Table  2.  Initial  values  and  perturbation  parameters. 


Penetrator 

Target 

L„ 

(m) 

Vi 

(m/s) 

e 

Tungsten 

Steel  (RHA) 

0.500 

2000 

0.0140 

Tungsten 

Steel  (RHA) 

0.500 

3000 

0.0063 

93%  tungsten 

Steel  (RHA) 

0.500 

2000 

0.0511 

93%  tungsten 

Steel  (RHA) 

0.500 

3000 

0.0227 

Depleted  uranium  (DU) 

Steel  (RHA) 

0.500 

2000 

0.0430 

Depleted  uranium  (DU) 

Steel  (RHA) 

0.500 

3000 

0.0191 

Forrestal  steel 

Silica-sand 

0.039 

3000 

0.0192 
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Figure  19.  Comparison  between  rod  and  penetration  velocities  for  a  first-order  perturbation 
solution,  a  third-order  perturbation  solution,  and  the  exact  Tate  solution.  This 
case  is  for  a  depleted  uranium  rod  impacting  a  steel  target  at  an  impact  velocity  of 
2000  m/s.  The  initial  rod  length  was  0.500  m. 


24 


Figure  20.  Penetration  depth  for  a  first-order  perturbation  solution,  a  third-order  perturbation 
solution,  and  the  exact  Tate  solution.  This  case  is  for  a  depleted  uranium  rod 
impacting  a  steel  target  at  an  impact  velocity  of  2000  m/s.  The  initial  rod  length 
was  0.500  m. 
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Figure  21.  Rod  length  for  a  first-order  perturbation  solution,  a  third-order  perturbation  solution, 
and  the  exact  Tate  solution.  This  case  is  for  a  depleted  uranium  rod  impacting  a  steel 
target  at  an  impact  velocity  of  2000  m/s.  The  initial  rod  length  was  0.500  m. 
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Figure  22.  Comparison  between  rod  and  penetration  velocities  for  a  first-order  perturbation 

solution,  a  third-order  perturbation  solution,  and  the  exact  Tate  solution.  This  case  is 
for  a  depleted  uranium  rod  impacting  a  steel  target  at  an  impact  velocity  of  3000  m/s. 
The  initial  rod  length  was  0.500  m. 
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Figure  23.  Penetration  depth  for  a  first-order  perturbation  solution,  a  third-order  perturbation 
solution,  and  the  exact  Tate  solution.  This  case  is  for  a  depleted  uranium  rod 
impacting  a  steel  target  at  an  impact  velocity  of  3000  m/s.  The  initial  rod  length  was 
0.500  m. 
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Figure  24.  Rod  length  for  a  first-order  perturbation  solution,  a  third-order  perturbation 
solution,  and  the  exact  Tate  solution.  This  case  is  for  a  depleted  uranium  rod 
impacting  a  steel  target  at  an  impact  velocity  of  3000  m/s.  The  initial  rod 
length  was  0.500  m. 


3.  Conclusions 


A  perturbation  solution  of  the  Alekseevski-Tate  equations  was  obtained  through  the  third  order. 
Agreement  with  the  exact  solution  is  excellent  for  a  tungsten  rod  impacting  a  steel  target  at  a 
velocity  of  2.0  km/s,  as  well  as  for  other  cases  of  interest  to  ballisticians  involved  in  the 
penetration  of  steel  targets.  Also,  the  velocities  and  residual  length  of  a  tungsten  rod  perforating 
a  finite-thickness  steel  plate  show  excellent  agreement  with  the  exact  equations.  The 
approximate  (perturbation)  solution  yields  the  dependent  variables  (V,  U,  X,  and  P)  as  explicit 
functions  of  time.  Also,  alternate  forms  of  the  normalization  of  the  pertinent  equations  are 
investigated  to  obtain  a  perturbation  parameter  much  less  than  one  for  various  penetration 
problems  depending  on  the  input  conditions,  namely  target  and  penetrator  densities,  strengths, 
and  initial  impact  conditions.  Also,  comments  are  made  regarding  the  singularity  in  time. 
Extensions  of  the  method  to  higher  orders  (i.e.,  the  fourth  order)  are  possible,  but  a  point  of 
diminishing  returns  has  been  reached.  The  current  third-order  solution  is  expressible  as  an 
algebraic  equation,  amenable  to  a  spread  sheet  or  simple  calculator  evaluation. 
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Appendix.  FORTRAN  Computer  Program 


Pertubation.f90 

FUNCTIONS: 

Pertubation  -  Closed-form  equations  for  Tate's  model  from  a  pertubation  analysis. 


PROGRAM:  Pertubation 

PURPOSE:  This  program  calculates  the  rod  velocity,  penetration  velocity,  penetration 

depth,  and  instantaneous  length  using  pertubation  method  up  to  the  third-order 
correction  terms. 

program  Pertubation 
implicit  none 
!  Indentifiers 


RT: 

Strength  resistance  of  target  (Kg/m.sA2) 

RHO  P: 

Density  ofpenetrator  (Kg/mA3) 

RHO  T 

Density  of  target  (Kg/mA3) 

YP: 

Flow  stress  ofpenetrator  (Pa) 

VI: 

Impact  velocity  (m/s) 

VR: 

Rod  velocity  (m/s) 

UP: 

Penetration  velocity  (m/s) 

LO: 

Original  length  of  penetrator  (m) 

LI: 

Instantaneous  length  of  rod  (m) 

P: 

Penetration  Depth  (m) 

T: 

Time  (microsecond) 

VN: 

Normalized  rod  velocity 

UN: 

Normalized  penetration  velocity 

PD: 

Normalized  penetration  depth 

LAMDA: 

Normalized  length 

TAU: 

Normalized  time 

MHU: 

Dimensionless  variable 

BETA: 

Dimensionless  variable 

ALPHA: 

Dimensionless  variable 

EPSI: 

Dimensionless  variable 

A,B: 

Expressions  used  for  equation  compactness 

DOUBLEPRECISION  RT,  RHO_P,  RHO  T,  YP,  VI,  VR,  UP,  LO,  LI,  P,  PD,  T,  VN,  UN,  LAMDA,  TAU, 
MHU,  BETA,  ALPHA,  EPSI,  VO,  UO,  LAMDAO,  PDO,  VI,  Ul,  LAMDA  1,  PD1,  V2,  U2,  LAMDA2,  PD2,  V,  U, 
CO,  Cl,  C2,  C3,  V3,  U3,  PD3,  LAMDA3,  THK,  A,  B 

INTEGER::  ORDER 

!  Input  values  for  velocity,  length,  and  material  constants 


31 


PRINT*,  "Enter  impact  velocity  of  penetrator,  'Meters  per  Second'" 

READ*,  VI 

PRINT*,  "Enter  original  length  of  penetrator,  'Meters'" 

READ*,  LO 

PRINT*,  "Enter  density  of  penetrator  material,  'Kilograms  per  Meter  Cubed'" 

READ*,  RHOP 

PRINT*,  "Enter  density  of  target  material,  'Kilograms  per  Meter  Cubed'" 

READ*,  RHOT 

PRINT*,  "Enter  flow  stress  of  penetrator  material,  'Pascals'" 

READ*,  YP 

PRINT*,  "Strength  resistance  of  target,  'Pascals'" 

READ*,  RT 

!  Check  for  data  inconsistencies 

!  The  "strength  resistance  of  the  target  (RT),"  cannot  be  equal  to  the  "flow  stress  of  the 

!  penetrator  (YP),"  for  the  pertubation  solution  to  be  meaningful. 

IF(RT.EQ.YP)THEN 

PRINT*,  "*****rj  must  not  equal  YP  for  a  meaningful  pertubation  solution.*****" 

STOP 

ENDIF 

PRINT*,  "Thickness  of  the  target,  'Meters'" 

READ*,  THK 

!  Assume  a  semi-infinite  plate  for  zero  thickness. 

IF(THK.EQ.O)THEN 

THK=1 

ENDIF 

PRINT*,  "What  order  pertubation  solution  do  you  require" 

READ*,  ORDER 

OPEN(UNIT=12,  FILE="Pertubation.dat",  STATUS="OLD") 

!  Normalization  scheme 

MHU=SQRT  (RHOT  /RHOP ) 

BETA=MHU/(  1 . 0+MHU)*  VI/LO 
ALPHA=(RT  -  YP)/YP 
EPSI=YP/(RHO_P*VI**2) 
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A=(1.0+MHU)/MHU 
B=  1 .0- ALPHA/MHU 


T=0.0 


DO 


TAU=BETA*T 

Zero  order  calculation  resulting  from  pertubation  analysis 
V0=1.0 

U0=1 .0/(1 .0+MHU) 

LAMDA0=1.0-TAU 
PD0=T  AU/MHU 

First  order  calculation  resulting  from  pertubation  analysis 


V1=((1.0+MHU)/MHU)*LOG(1.0-TAU) 

Ul=(  1 ,0/MHU)*(LOG(  1 .0-TAU)-ALPHA) 

LAMDA1=((1.0+MHU)/MHU)*((1.0-ALPHA/MHU)*TAU+(1.0-TAU)*LOG(1.0-TAU)) 

PD1=-(1.0+MHU)/(MHU**2)*(TAU*(ALPHA+1.0)+(1.0-TAU)*LOG(1.0-TAU)) 


Second  order  calculation  resulting  from  pertubation  analysis 

V2=((  1 ,0+MHU)/MHU)**2*(  1 ,0-ALPHA/MHU)*(LOG(  1 .0-TAU)+l/(l  .0-TAU)- 1 .0)-  & 
0.5*(((1.0+MHU)/MHU)**2)*(LOG(1.0-TAU))**2 

U2=V2/(  1 ,0+MHU)+(ALPHA/MHU**2)*(  1 ,0+MHU)*LOG(  1 ,0-TAU)+(ALPHA**2)*(  1 .0-& 

MHU*  *2)/(2. 0*MHU*  *3 ) 

C0=3.0-3.0*ALPHA/MHU+ALPHA**2/(2.0*MHU**2)-ALPHA**2/(2.0*MHU) 

LAMDA2=((1.0+MHU)/MHU)**2*(-(3.0-3.0*ALPHA/MHU+ALPHA**2/(2.0*MHU**2)-& 
ALPHA**2/(2.0*MHU))*(l-TAU)+(1.0-ALPHA/MHU)*LOG(1.0-TAU)+2.0*(1.0-& 
ALPHA/MHU)*(1.0-TAU)*LOG(1.0-TAU)-0.5*(1. 0-TAU)*  (LOG(1.0-TAU))**2+C0) 

PD2=(1.0+MHU)**2/MHU**3*(B*(1.0-2.0*TAU-(1.0-TAU)*LOG(1.0-TAU)-LOG(1.0-&  TAU))+(1.0- 
TAU)/2.0*(LOG(  1 .0-TAU))**2-(  1 ,0-TAU)*LOG(1 .0-TAU)+ 1 .0-&  TAU+ALPHA*(  1 .0-TAU-(l  .0- 
TAU)*LOG(  1 .0-TAU))+  ALPHA**2*(1.0-&  MHU)*TAU/(2.0*MHU)-(2+ALPHA-ALPHA/MHU)) 


Third  order  calculation  resulting  from  pertubation  analysis 
Cl=-3.0/2.0*B**2-2.0*B+ALPHA**2/(2.0*MHU)*(1.0-1.0/MHU) 

V3=A**3*((4.0-5.0*ALPHA/MHU-ALPHA**2/(2.0*MHU)+3.0*ALPHA**2/(2.0*MHU**2))& 

*LOG(1.0-TAU)-0.5*B**2*1.0/(1.0-TAU)**2+(2.0*B+2.0*B**2-ALPHA**2/(2.0*MHU)+& 
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ALPHA**2/(2.0*MHU**2))*1.0/(1.0-TAU)-2.0*B*(LOG(1.0-TAU))**2+0.5*(LOG(1.0-&  TAU))**3- 
B*(LOG(  1 . 0-TAU))/(  1 .0-TAU)+C  1 ) 

U3=(V3*(V0-U0)+(V1-U1)*(V2-U2)-MHU**2*U1*U2)/((V0-U0)+U0*MHU**2) 

C2=(1.0+MHU)**3/MHU**4*(-4.0-7.0/2.0*B**2-7.0*B+9.0/2.0*ALPHA**2/MHU-&  5.0*ALPHA- 
3.0*ALPHA**2/2.0-ALPHA**3/(2.0*MHU**2)*(1.0-MHU)**2+B**2/2.0) 

PD3=-(1.0+MHU)**3/MHU**4*((1.0-TAU)*(-4.0-7.0/2.0*B**2-7.0*B+9.0/2.0*&  ALPHA**2/MHU- 
5.0*ALPHA-3.0*ALPHA**2/2.0-ALPHA**3*& 

(1.0-HU)**2/(2.0*MHU**2))  +(1.0-TAU)*LOG(1.0-TAU)*(7.0/2.0+6.0*B+& 
+3.0*ALPHA*B+3.0*B**2/2.0+ALPHA+3.0*ALPHA**2/2.0)+B**2/(2.0*(1.0-TAU))  & 
+(B+5.0/2.0*B**2+ALPHA-3.0*ALPHA**2/(2.0*MHU)+1.0/2.0)*LOG(1.0-TAU)+(1.0-& 
TAU)*(LOG(  1 .0-TAU))**2*(-2.0*B-3 .0/2.0*(  1 .0+ALPHA))-B/2.0*(LOG(  1 .0-TAU))**2+(  1 .0-& 
TAU)/2.0*(LOG(1.0-TAU))**3)+C2 

C3=((l+MHU)/MHU)**3*(5.0/2.0+13.0/2.0*B**2+5.0*B-5.0*ALPHA**2/(2.0*MHU)-& 

ALPHA**3/(2.0*MHU**3)*(MHU-1)**2) 

LAMDA3=((1+MHU)/MHU)**3*((-7.0*B/2.0)*(1.0-TAU)*(LOG(1.0-TAU))**2+(7.0*&  B**2/2.0  - 
ALPHA**2/(2.0*MHU)+1.0/2.0)*LOG(1.0-TAU)  -B/2.0*(LOG(1.0-TAU))**2&  +(2.0+5.0*B+4*B**2- 
2*ALPHA**2/MHU)*(1.0-TAU)*LOG(1.0-TAU)+(-5.0/2.0-7.0*B**2& 
-5.0*B+5.0*ALPHA**2/(2.0*MHU)+ALPHA**3/(2.0*MHU**3)*(MHU-1)**2)& 
*(1.0-TAU)+B**2/(2.0*(1.0-TAU))+(1.0-TAU)/2.0*(LOG(1.0-TAU))**3)+C3 

Summation  of  series  for  1st  order,  2nd  order,  and  3rd  order  pertubation  solution 


IF(ORDER.EQ.l)THEN 

V=V0+EPSI*V1 

U=U0+EPSI*U1 

LAMDA=LAMDAO+EPSI*LAMDA  1 
PD=PDO+EPSI*PD  1 
ELSEIF(ORDER.EQ.2)THEN 

V=V0+EPSI*V1+(EPSI**2)*V2 

U=U0+EPSI*U1+(EPSI**2)*U2 

LAMDA=LAMDA0+EPSI*LAMDA1+(EPSI**2)*LAMDA2 

PD=PD0+EPSI*PD1+(EPSI**2)*PD2 

ELSEIF(ORDER.EQ.3)THEN 

V=V0+EPSI*V1+(EPSI**2)*V2+(EPSI**3)*V3 

U=U0+EPSI*U1+(EPSI**2)*U2+(EPSI**3)*U3 

LAMDA=LAMDA0+EPSI*LAMDA1+(EPSI**2)*LAMDA2+(EPSI**3)*LAMDA3 
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PD=PD0+EPSI*PD1+(EPSI**2)*PD2+(EPSI**3)*PD3 

ENDIF 

!  Converting  normalized  values  into  actual  values 

VR=V*VI 

up=u*vi 

P=PD*LO 

LI=LAMDA*LO 

!  Penetration  process  stop  conditions 

!  For  RT  >  YP,  the  minimum  rod  velocity  required  to  achieve  penetration  occurs  when  the  penetration 
!  velocity  is  zero  (UP  =  0).  At  this  point  the  target  begins  to  behave  as  a  rigid  body  and  the  rod 
!  bounces  off  the  target. 

IF(RT.GT.YP.AND.VR.LE.SQRT(2.0*(RT-YP)/RHO_P))EXIT 

!  For  RT  <  YP,  the  minimum  rod  velocity  required  to  achieve  penetration  occurs  when  the  penetration 
!  velocity  is  equal  to  the  rod  velocity  (UP  =  VR).  At  this  point  the  rod  begins  to  behave  as  a 
!  rigid  body. 


IF(RT.LT.YP.AND.VR.LE.SQRT(2.0*(YP-RT)/RHO_T))EXIT 

!  If  the  penetration  is  greater  than  the  thickness  of  the  target,  perforation  is  achieved  terminate 
!  loop 


IF(P.GT.THK)EXIT 

!  If  the  instantaneous  length  is  less  than  zero,  then  the  rod  had  totally  eroded,  Terminate  Loop 
IF(LI.LT.0.0)EXIT 
!  Output  data  for  graphical  processing 

WRITE(12,  100),  T,  VR,  UP,  P,  LI 

100  FORMAT(F10.7,  5X,  F15.10,  5X,  F15.10,  5X,  F10.8,  5X,  F10.8) 

!  Time  increment,  1  microsecond  step  size 
T=T+0. 000001 
END  DO 
PRINT*,  EPSI 

END  PROGRAM  Pertubation 
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